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S : ABSTRACT 

Using the most recent results about white dwarfs in 10 open clusters, we 
revisit semi-empirical estimates of the initial-final mass relation in star clusters, 
with emphasis on the use of stellar evolution models. We discuss the influence 
of these models on each step of the derivation. One intention of our work is to 
use consistent sets of calculations both for the isochrones and the white dwarf 
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cooling tracks. The second one is to derive the range of systematic errors arising 
from stellar evolution theory. This is achieved by using different sources for the 
stellar models and by varying physical assumptions and input data. We find that 
systematic errors, including the determination of the cluster age, are dominating 
the initial mass values, while observational uncertainties influence the final mass 
primarily. After having determined the systematic errors, the initial-final mass 
relation allows us finally to draw conclusions about the physics of the stellar 
models, in particular about convective overshooting. 

Subject headings: stars: AGB and post-AGB — stars: evolution — stars: indi- 
vidual(Sirius) — stars: mass loss — white dwarfs — open clusters and associa- 
tions: individual (Pleiades, Hyades, Praesepe, NGC1039, NGC2168, NGC2099, 
NGC2516, NGC3532, NGC6819, NGC7789) 



1. Introduction 

The initial-final mass relation (IFMR) for low- and intermediate-mass stars is an im- 
portant input for many astrophysical problems. Given the initial Main Sequence (MS) mass 
of a formed star, the IFMR provides the expected mass during its final White Dwarf (WD) 
cooling stage, and is an estimate of the total mass lost by the star during its evolutionary 
history. A correct assessment of the IFMR is very important when predicting, for example, 
the chemical evolution history of stellar populations, or their mass-to-light ratio (defined as 
the ratio of the mass of evolving stars plus remnants - WDs, neutron stars and black holes 
- to the integrated luminosity of the population), and in general for any problem related 
to the origin and evolution of gas in stellar populations. It is also a crucial item when 
using the WD luminosity functi ons to determine the age of stellar populations (see, e.g.. 



Prada Moroni fc Stranierd 120071 ). It also provides an empirical test concerning the upper 
initial mass value Mup of stars developing degenerate carbon-oxygen cores, and thus a lower 
limit of the mass range of core-collapse super novae. For a thorough review of the history of 



the IFMR, we recommend IWeidemannI (J2000l ) 



Theoretical estimates of the IFMR are still prone to large uncertainties. This is due 
to our poor knowledge of the efficiency of mass loss processes for low- and intermediate- 
mass stars but also to uncertainties in the predicted size of CO cores during the Asymptotic 
Giant Branch (AGB) evolution, resulting mainly from uncertainties in the treatment of the 
thermal pulse phase, with the associated third dredge- u p, hot bottom burning (see, e.g.. 



Iben &: Renzinilll983l ) and also the treatment of rotation (JDominguez et al.lll996l ). In fact 



the recognition that intermediate-mass stars may lose the largest part of their mass during 
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the AGB phase stems from early IFMRs ( Koester &: Weidemannlll980l ) 



Starting with the pioneering work by IWeidemannl (119771 ). (semi-)empirical routes have 
been followed to establish the IFMR independently of theoretical mo delling the AGB phase, 
l eadiri g to a series o f global determinations o f the IFMR by, e.g.. IWeidemann fc Koester 
(1l983l ). IWeidemanru (119871 ). IWeidemannI (120001 ). iFerrario et al.l (120051 ). and many others. 



Methods to estimate the IFMR that make use of the smallest number of assump- 
tions are probably t hose based on the use of WDs harbored in star clusters (but see also 
Catalan et al.ll2008bl . for a very recent study of the IFMR based on WDs in common proper 
motion pairs) . Thanks to a large amount of observational effort by various authors (see, e.g . 



Koester &: Reimers 1993 



Barstow et al. 2005: iDobbie et al 



1996 



C 



aver et al. 2001: Dobbie et al. 2004: Kahrai et al 



2006al ) the recent study of the IFMR by iFerrario et al. 



2005 



( 120051 ) makes use of 40 DA ( hydro gen-atmosphere) WDs belonging to 7 open clusters. Even 
more recently iKalirai et al.l (120081 ) have ad ded a few more data points below the low-mass 
end (initial masses below ~2.0 Mq) of the IFerrario et al.l (120051 ) sample, by in cluding WDs 



detect ed in the old open clusters NGC6819, NGC7789 and NGC6791, while iRubin et al. 
J2008h have provided WD data for the ~200 Myr old cluster NGC1039. Therefore, the ob 



servational material is increasingly covering not only more objects, but also a larger range of 
cluster - and thus WD - ages and metallicities, such that even differential investigations may 
become feasible in the near future. At the same time, high-resolution, high signal-to-noise 
spectra and improved spectral analysis methods reduce the intrinsic errors in the determi- 
nation of the WD masses. It is therefore time to investigate the systematic errors arising 
from stellar evolution theory and models that enter these semi-empirical methods. This is 
the main object of this work. 

IFMR determinations based on cluster WDs work as follows: after detection, spectro- 
scopic estimates of the WD surface gravity g and Teg are needed. For a fixed g-T^s pair, 
interpolation within a grid of theoretical WD models covering a range of masses provides the 
mass Mf and the age tcooi of the WD. Independent theoretical isochrone fits to the Turn Off 
(TO) luminosity in the cluster Color-Magnitude-Diagram (CMD) provide an estimate of the 
cluster age tcius- The difference tcius — ^cooi is equal to the lifetime tprog of the WD progenitor, 
from the MS until the start of the WD cooling. Making use of mass-lifetime relationships 
from theoretical stellar evolution models, the initial progenitor mass Mi is immediately ob- 
tained from tprog. 

It is clear how and where theoretical models play an important role in this procedure. 
WD cooling models plus progenitor evolutionary tracks and isochrones are needed to deter- 
mine the various tdus, tcooi and tprog- Given the very short timescale of the AGB evolution 
(see also next section for a short discussion on this issue) AGB modelling and its related 



-4- 



uncertainties do not play a role in these types of analyses, but uncertainites entering the 
previous evolutionary phases and the WD cooling stage can affect the derived IFMR. An- 
other important point is the issue of consistency. First of all, the isochrones employed to 
determine tcius should be computed starting from the same evolutionary tracks adopted to 
determine Mi from the estimated value of tprog- Consistency is not usually achieved in ex- 
isting estimates of the IFMR. Even in the most recent studies of the IFMR, in several cases 
cluster ages are obtained from earlier studies that have employed models different from the 
ones used to estimate M^. Secondly, all cluster ages should be determined following the same 
method, and - most importantly - employing the same set of isochrones. Instead, results 
from different sources for the various clusters are used. 



As an example, we s ummarize briefly the w ork by iFerrario et al.l (120051 ). The Hyades 
age has been taken from IPerryman et al.l ( 19981), who had calculated their own isochrones 
with an upda ted version of the CESAM code (JMorell 119971 ). In case of M35 the age quoted 
in Table 1 of Ivon Hippell (120051 ) was employed; this age was derived summarizing the re- 
sults of three different previous works by other authors. The same is true for the Pleiades, 
where the assumed cluster age comes from a synthesis of the results of sever al groups, ob 



tained wit h different generation s of stellar evolution models. The isochrones by lSchaller et al. 



(Il992h led iMevnet et al.l (119931 ) to an age of 141 Myrs for NGC2516, which was quoted by 
Ferrario et al.l (120051 ) as 158 Myr. The remaining clusters in that sample also have ages 
from various isochrone sources, which span more than a decade of progress in stellar evolu- 
tion modelling, most notably the sig nificant improv e ments in opacities and the equation of 
state. The cooling t racks came from iFontaine et al.l (120011 ) . and the tprog^ M^ relation from 
Girardi et al.l (120021 ). taken from the tracks with the metallicity closest to the adopted one, 
but without interpolation between compositions. It is obvious that there is a certain degree 
of arbitrariness in the use of stellar evolution theory. This inherent inconsistency is probably 
the result of the fact that most works on the IFMR are more focused on acquiring new and 
improved data with an emphasis on the internal errors of their derived T^s and Mf. 

The aim of the present work is instead to concentrate on the stellar evolution models 
needed for deriving the IFMR. It is obvious from the previous discussion that a consistent 
use of stellar models and procedures about how to employ them is required, to achieve a 
higher accuracy. We therefore will make use of calculations of stellar models, that cover 
all phases from the MS to the WD cooling phase. They will be used for the appropriate 
metallicity of each cluster, in order to have at hand accurate and homogeneous isochrones 
and progenitor lifetimes. This will provide an IFMR with the highest internal consistency 
achieved so far. 



To get a measure of the errors introduced by using different codes and input physics. 
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we also repeat the procedure by using alternative sets of models for the WD cooling and/or 
progenitor and cluster ages. While this gives a first indication about the size of the systematic 
errors resulting form different stellar model sources, we finally quantify the influence of 
various physical effects important for both the isochrone and WD cooling calculations. Using 
Monte Carlo methods we determine the total extent of the systematic errors. This allows 
us also to verify which features of the IFMR are robust and which ones depend on present 
model uncertainties. 

Using the robust results, we can finally exclude some of the stellar models because they 
lead to internal contradictions. One of our conclusions is that models without overshooting 
from convective cores on the MS do not yield self-consistent results. On the other hand, the 
parametrization of overshooting during the AGB phase, which influences the growth of the 
carbon-oxygen core, cannot be the same in all convective layers as it is at the convective 
core during the MS phase, because it produces final masses inconsistent with observations. 
We thus can also reduce uncertainties in the physics of stellar interiors. 

The structure of the paper is as follows: In Sect. 2 we outline the procedure to obtain 
an IFMR, summarize the observational input data we have used for the WDs and the parent 
clusters, and present the set of the stellar models used for determining the reference IFMR 
as well as for the evaluation of the systematics. Section 3 then concentrates on the estimates 
of the global error bars, while Sect. 4 contains the results and a comparison with theoretical 
IFMRs. A summary and conclusions follow in Sect. 5. 



2. IFMR determination 

As discussed in the Introduction, the method we follow to determine the IFMR makes 
use of WDs in star clusters. Observations have to provide us with: 

• Surface gravity g and effective temperature T^^ of individual WDs. 

• Observed CMDs, estimates of [Fe/H] and E(B — V) for the clusters harboring the 
observed WDs. 

Theory has to provide: 

• Sets of stellar evolutionary tracks, plus isochrones - transformed to the CMD - for a 
range of age and [Fe/H] values, to determine - in conjuction with theoretical, empir- 
ical or semi-empirical distance determination methods - the cluster ages, plus stellar 
lifetimes. 
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• WD cooling tracks covering a range of masses, to provide the cooling ages of the 
observed WDs. 

Concerning the cooling tracks, the metallicity of the WD models is irrelevant as the 
envelopes are of pure H (He and metals have settled below the envelope due to diffusion) 
and the previous core composition may influence only the chemical profile of the core. This 
will be investigated in a broader context later in this study. 

The method works as already summarized in the Introduction; for each cluster WD 
with an estimated pair of g and Tgfj values, linear interpolation within a grid of theoretical 
WD models provides the mass Mf (mainly fixed by the value of g) and the age tcooi (mostly 
determined by the value of Teg) of the WD. Isochrone fits to the TO region of the observed 
V — (B — V) CMD determine the cluster age tcius- If the actual cluster [Fe/H] does not 
correspond to any of the values in the isochrone grid used, we determined the age with the 
two closest grid values bracketing the cluster [Fe/H]. The final cluster age is then obtained 
by linear interpolation between these two values. 

The difference tcius — ^cooi represents the lifetime tprog of the WD progenitor, from the 
MS until the start of the WD cooling. Theoretical mass-lifetime relationships for the cluster 
metallicity, obtained from the same evolutionary calculations, finally provide the initial pro- 
genitor mass Mi. We employ a quadratic interpolation in mass and a linear interpolation in 
[Fe/H] among the available model grids to determine Mi from the value of tprog- 

A potential inconsistency in this procedure is the effect of the duration of the AGB ther- 
mal pulse phase on the estimate of M^. Theoretical models without inclusion of the thermal 
pulse phase, or predicting an IFMR different from the one obtained semi-empirically, might 
be providing inconsistent thermal pulse lifetimes, hence wrong values of the total lifetimes, 
that in turn affect the derivation of M^. Luckily, this effect is negligible; estimates of ther 



mal pulse lifetimes obtained f rom synthetic AGB modelling (jWagenhuber fc Groenewegen 



19981 : iMarigo fc Girardil 120071 ) show that the duration of this phase is typically at most 
within a few Myr (decreasing sharply with increasing total mass) even in extreme cases of 
the largest possible growth of the CO core. In fact, from full evolutionary stellar models. 



Serenelli &: Fukugital (120071 ) have shown that the duration of the AGB phase (including both 
the early- AGB and the thermally pulsating phases), for stars with initial masses between 1 
and 8 M0, is at most 1% of the MS lifetime. 
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2.1. WD data 



Our reference data set comprises 52 WDs in 10 open clusters, plus Sirius B, for a 
total of 53 DA objects, as re ported in Tab l e [H We considered the same cluster sam- 
ple as in the recent work b y iFerrario et al.l (120051 ) , that contains data for WDs in th e 



Pleiades (IFerrario et al.ll2005[). Hvad es and Praes epe (Claver et al.ll200ll: iDobbie et al.ll2004). 



20051), NGC2168 (IWilliams et al. 



2004; Ferrario et al. 



200 



K 



NGC 2516 JKoester fc Reimerslll996h . N GC3532 jKoester fc ReimerJ l993h. NGC2099 JXahrai et al 



We wish to notice that our 



Pleiades sample contains only the object listed by IFerrario et al.l (120051). We did no t includ e 
data for additional two WDs (GD 50 and PG 0136+251) discussed bv bobbie et all J2006bh . 
because their membership of the Pleiades is not absolutely certain. 

To this cluster sample we added NGC6819, NGC7789 and NGC10 3 9. Fo r the WDs 
in the first t wo cl usters, Tes and g values are taken from iKalirai et al.l (120081 ) . and from 
Rubin et al.l (120081 ) for the last one. 



The TefT, g and corresponding error estimates for the clusters in common with lFerrario et al 



( I2OO5I ) are the same as in their Table 1, with a few exceptions. Values of Tes and g for two 
Praesepe WDs (0837+ 185 and 0837+218), id entified as WD #13 and #14 in Table H have 
been updated following iDobbie et al.l (l2006al ). The same authors provide data for four addi- 
tional WDs in the same cluster (ID numbers from 16 to 19 in Table H]), that we have included 
in our analysis. And finally, we also use their data for the single Pleiades WD. 

As an additional note about the cluster WD sample, our referee (K. Williams) has 
pointed out that the Hyades WD 0437 +138, that appears as WD #8 in Table [U has been 



incorrectly included as a DA object in (IFerrario et al.ll2005l ) study. In fact, it is a DBA WD, 
with log(n///nj:/e)=— 4.5 in the atmosphere. To assess the errors in the derived mass and 
cooling age for this object when using DA models, we have made the following test. We 
have compared log{g) values derived from DB and DA models at the observed Tefj of WD 
#8, for the two extreme cases of 0.55 and 1.0 Mq (the mass we obtain from the analysis 
with DA models is ~O.77M0). Our adopted DB models have been computed as described 



m 



Salaris et all (I2OOII ). 



For both masses the difference in log{g) derived from the DA and DB tracks is smaller 
than the la error bar on the observed value for WD #8 (as listed in Table [T]). Also the 
difference in cooling ages at the observed Tgfj is well below the error bars we derive from 
our analysis, as reported in Table [51 This implies that the analysis of WD #8 with DB 
models would provide mass and cooling age compatible (within the la errors) with the 
results obtained from the DA models employed in this analysis. The reason is that the 
mass-radius relation is not strongly affected by the atmospheric composition, and also at 



the relatively high Tcs of this object the evolutionary timescales of DA and DB models are 
similar. Moreover, given that this WD does not have a pure-He envelope, differences will be 
even more negligible. 

In addition to a large sample of cl uster WDs, we also include Sirius B in our analysis. 
The data for Siri us B (also included inlFerrario et al.ll2005l ) in Table [1] are from the recent 
determination by lBarstow et al.l (120051 ). From the estimated mass and radius of the primary 
MS star (Sirius A) one can derive its age by comparison with stellar evolut ionary tracks 



i n the mass-radius plane. We assumed [Fe/H]=0.0 for Sirius A, following iLiebert et al. 



(I2OO5I ). to which we associated an 0.10 dex la error. The age of the star is obtained from 
models bracketing the metallicity of the system (if different from the values provided by the 
employed isochrone grid). Linear interpolation in [Fe/H] between the two bracketing results 
gives then the age for the appropriate [Fe/H]. The age of Sirus A provides the age of the 
system tgys- Sirius B cooling time tcooi and mass are determined from the spectroscopic Tes 
and g estimates, as described before. The difference tgys-tcooi provides the age of the Sirius B 
progenitor, hence Mi after using a grid of theoretical lifetimes for solar metallicity. 



2.2. Cluster data 



The sources for the cluster [Fe/H], E(B-V), and their CMDs used for age determinations, 
are listed in Table [2l 



Values for [Fe/H] have been taken mainly from the compilation by iGrattonI (J2000l ). who 
presents a collection of [Fe/H] estimates for a large sample of open clusters, recalibrated 
against results from high-resolution spectros c opy. Alternative estimates have been adopted 
for NGC2099, which is not in the ICrattonl J2000h compilation, NGC6819 and NGC1039, 
for which recent high-resolution spectroscopic measurements have appeared in the literature 
(see Table 12]). 



Reddenings are taken from the compilation by lLoktin et al.l (1200 ll ). These are weighted 
averages of estimates obtained with various methods, all based essentially on color-color 
relationships. For NGC6819 we em ployed a very recent E(B-V) determination based on 
spectroscopy (JBragaglia et al.l I2OOII). and for th e Plei ades reddenings to individual stars 
have been adopted, as discussed in IPercival et al.l (120031 ). The case of NGC7789 is discussed 
below. 

In case of Hyades and Praesepe we have assumed - as usual - zero reddening. Through- 
out this paper we use Ay = 3.2 E(B — V). We adopted a 0.02 mag uncertainty in E(B — V) 
when no error is quoted by the source of the data. 
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To estimate the age of the whole cluster sample we have followed a two-step procedure 
where we tried to minimize the inputs from theoretical isochrones. First, cluster distance 
moduli have been determined using a fully empirical MS-fitting to the V — (B — V) CMDs, 
that makes u se of the sampl e of fie ld dwarf with accurate Hipparcos parallaxes and [Fe/H] 
presented by iPercival et al.l (120031 ) . Details about the method can be found in the same 
paper. 

Here we just recall that for each cluster a MS fiducial line has been derived, by plotting 
color histograms in V-magnitude bins, and considering the mode of the resulting distribution. 
In this way one minimizes the impact of unresolved binaries. The resultant points are fitted 
to a cubic function that provides the cluster fiducial line. After applying (empirical) color 
corrections that account for the actual cluster metallicity and for the reddening E(B-V), the 
field dwarf template sequence is shifted in magnitude to match the MS fiducial line, and the 
apparent distance modulus is obtained by x^ minimization. For objects like NGC2168 or 
the Pleiades, with an extremely well defined and thin (in color) MS, the separation of the 
binary sequence is very evident, and a cubic fit directly to the single star MS is performed 
in order to derive the fiducial line. 

In case of NGC3532 we could not find available data with good photometry in the MS 
magnitude range of our template field dwarf sample. The distance modulus of this cluster 
has been obtained by comparison with Praesepe, that has the same [Fe/H] within errors 
(and also a comparable age). The larger (as compared to our other estimates) error bar on 
its distance (hence on its age) reflects this less direct determination of the distance modulus. 



As a final remark, we wish to notice that, as discussed i.e. in iPercival et al.l (120031 ). 
a simultaneous MS-fitting to both VI and BV data can provide both the cluster distance 
modulus and reddening. This procedure has been applied (using the same field dwarf sample, 
th at has homogeneo us BVI photometry) to NGC7789, employing additi onally the VI data 



bv iGim et al.l (jl998|). Our adopted E(B-V) for NGC7789, used also in IPercival fc Salaris 
( 120031 ) . has been determined in this way. This method to determine the reddening could not 
however be applied to the whole cluster sample because of the lack of suitable BV and VI 
photometry for all clusters. 

For a few clusters we have used distances already published in the hterature, obtained 
with the same technique and the same Hipparcos field dwa rf sample. More in detail, for 
the Hyades we have used the res ult bv IPercival et al.l (120031) . that is in perfect agreement 
with the Hipparcos distance by iPerryman et al.l (Il998l). For the Pleiades we have used 



the distance modulus determined bv IPercival et al.l (120051 ). while for NGC7789 we have 
employed the result by 
of the investigation by 



Percival fc SalarisI (l2003h . The cluster NGC6791, which is also part 



Kalirai et al.l (120081 ) was not considered in our analysis because its 
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metallicity [Fe/H]~+0.45 (JGratton et al.ll2006l ) is higher than the upper [Fe/H] range of the 
field dwarf sample ([Fe/H] =+0.30), and therefore we could not apply this empirical method 
to determine its distance. Since we did not want to introduce a second, inconsistent method, 
we disregarded this cluster. 

Second, isochrone fits to the TO region (in the V — (B — V) plane), employing the 
distance moduli derived in the previous step, determine the cluster age. 

The final ages for various choices of the isochrone grid, as well as the adopted cluster 
[Fe/H] and E(B — V) values are reported in Table [31 Age uncertainties are due to both 
distance (which is in turn affected by the assumed [Fe/H] and E(B — V) values and their 
errors) and [Fe/H] uncertainties. Figure [1] compares the ages obtained with the BaSTI 
overshooting isochrones (that throughout the paper will be considered as ou r reference set 



of mo dels to determine the IFMR; see Sect. 2.3) with the ages adopted by iFerrario et al. 



( 120051 ) for the clusters in common (plus Sirius A). In spite of the dishomogeneity in terms 
of methods, data and isochrones employed, the two sets of ages compare well, within the 
respective error bars. The only exception occurs for NGC2099 (M37), for which our derived 
distance modulus is 0.4 magnitud es larger than i mplicit ly assumed in iFerrario et al.l (120051 ) 
- who adopt the age estimate by iKalirai et al.l (l2001bl ) - leading to a smaller cluster age 
(320 Myr as cor npared to 520 Myr). However, if we adopt a lower metallicity and reddening 
as suggested by iKalirai et al.l (120051 ) and discussed in Sect. 2.4 (see Table pi) we derive an 
age of 550 Myr for NGC2099, in close agreement with IFerrario et al.l (120051 ). Obviously, at 
this time the uncertainties in the cluster parameters dominate the determination over those 
from using different theoretical isochrones, which constitutes our first result. 



2.3. Models 



Our reference IFMR is defined by the use of the following set of models. We refer the 
reader to the original papers quoted below for details of the models and calculations. 



Pietrinferni et al.l (120041 ) models from the MS to the AGB with the inclusion of core 



overshooting during the central H-burning phase (hereafter BaSTI models with overshoot- 
ing), to estimate in a consistent way both t^ius and tprog- 



Salaris et al.l (120001 ) WD cooling models (hereafter SOO WD models) to estimate tcooi- 



The code employed to compute these WD models is the same as for the progenitors. Some 
major sources of input physics in the WD models (see SOO) are different from what has been 
employed in BaSTI models, most notably equation of state and low-temperature opacities. 
This is somewhat unavoidable, given the different range covered by many physical parameters 
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throughout WD structures, as compared to their progenitors. It is very important to remark 
here about the internal CO stratification of the WD models. The CO profile in the inner core 
of a WD is built up during the central He-burning phase, whereas the more external part 
of the WD core is processed during the early AGB phase and the thermal pulses. Current 
uncertainties in the treatment of mixing in stellar interiors and in some key reaction rate 
(i.e., the C^^ + a reaction) introduce additional uncertainties in the predicted CO profiles. 
SOO models employ a CO stratification obtained from evolutionary model s of solar metallicit y 



progenitors, from the MS until the completion of the first thermal pulse (jSalaris et al.lll997l ). 
The core mass and associated CO profile at the end of the pulse were taken as representative 
of the final WD object. When estimating the total error budget in the IFMR, we have 
accounted for the uncertainties in the CO profiles by using WD models with different C/0 
mass fraction ratios in their cores. This will be described in detail in Sect. [31 We anticipate, 
however, that these uncertainties are very small since all the WD stars in our sample have 
not entered yet the crystallization phase, where the detailed core composition could have a 
larger effect. 

We have studied systematic effects on the IFMR relation due to: 



i) Different model grids. We redetermine the IFMR employing iGirardi et al.l (120001 ) 
models from the MS to the AGB (hereafter Padua models; they also include core overshooting 
but with a different formalism compared to BaSTI), keeping the WD models unchanged 
(SOO). 



ii) Effect of overshooting. We redetermine the IFMR employing iPietrinferni et al.l (120041 ) 
models from the MS to the AGB without core overshooting (hereafter BaSTI no-OV models), 
again keeping the WD models unchanged (SOO). 

iii) Different WD model grids. We redetermine the IFMR employing the WD models 
co mputed with the LP CODE (hereafter LPCODE models; details of the LPCODE are given 
in lAlthaus et al.ll2003l ). keeping the progenitor models unchanged (BaSTI with overshooting). 
This tests not only different codes, but also - in an unsystematic way - variations of the 
physics employed and of the WD chemical stratification. 

iv) Systematic variation in WD input physics and chemical profiles. We redetermine 
the IFMR employing the LPCODE WD models, by varying the CO-core stratification, H- 
envelope thickness, electron conduction opacities and neutrino energy loss rates. In this case, 
the progenitor models are unchanged throughout (BaSTI with overshooting). 
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2.4. The reference IFMR 



We present in this section our reference IFMR, while we defer to subsequent sections 
the discussion of the uncertainties involved in the IFMR determination. We use the SOO WD 
models to derive the WD (final) masses and cooling ages. Using our reference cluster ages 
(fifth column in Table [3]) and the BaSTI models with overshooting we construct the IFMR 
by deriving the progenitor (initial) masses. Figure [2] displays this reference semi-empirical 
IFMR. Error bars include all uncertainty sources we consider. We defer a detailed discussion 
to Sections [3] and HJ but note here that the most relevant sources of uncertainty are due to 
the following effects: (i) errors in the WD T^s and g values; (ii) uncertainties in the cluster 
[Fe/H]; these affect the age through their infiuence on the derived cluster distance (mainly) 
and on the isochrones to be used for the age determination; (iii) errors in the reddening 
E(B — V), which infiuence the age determination through the effect on the cluster distance. 
In case of one star (WD #38 in NGC2099) the progenitor age results to be negative. We 
did not consider this object in the analysis of the IFMR. 

We have performed at this stage a simple check of ou r reference IFMR, comparing the 



distribution of Mj and Mf values with iFerrario et al.l (120051 ) - disregra ding the attached error 



bars - by means of a bi-dimensional Kolmogorov-Smirnov (KS) test (jPress et al.lll992l ). We 
considered only clusters (plus Sirius) in common between the two studies. As customary, we 
accept the existence of a significant difference between the two IFMRs if the probability P 
that the two samples of (Mi, Mf) pairs were drawn from different distributions was higher 
than 95%. We obtain P values smaller than 50%, highlighting the statistical agreement 



between our reference IFMR and lFerrario et al.l (120051 ) results for the clusters (plus Sirius B) 
in common. 

We have also performed analytical fits to the data. In doing this, we have left out 
the data from NGC2099 because our results for this cluster show some anomalies, e.g. a 
progenitor with negative age (WD 7,^38 mentioned above), a 0.45 Mq WD with a progenitor 
mass above 8 M© (WD #37) and some other issues to be discussed later in the paper. 

In Figure [2] we show two analytic fits to the data. The first one, linear, gives 

Mf = 0.084Mi + 0.466 (1) 

and a reduced chi-squared Xr = 3.7. Approximate values for the la rms dispersion about 
this relation are: 0.075 M© below 2 Mq, 0.12 Mq between 2.7 and 4 M©, 0.09 M© between 
4 and 6 Mq and then it increases almost linearly up to 0.20 Mq at around 8.5 Mq. 

Guided by theoretical models, which show a break in the slope of the IFMR at around 
4 Mq (see Fig. [3]), we have also performed a piecewise linear fit with a pivot point at 
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Mi = 4 M0 that yields 



0.134Mi + 0.331 1.7 M© < M; < 4 M© 



^ , w.x.,..., -r W...X X., .VX0 ^ ..., ^ , .VX0 

^ \ 0.047Mi + 0.679 4 M© < M;. ^ ^ 

This fit gives Xr = 2.7 at the expense of introducing only one more free parameter (we force 
the fit to be continuous at the pivot point). In this case, the dispersion is reduced to 0.05 M© 
below 2 Mq and to 0.09 Mq above 6 Mq, while it remains basically unchanged anywhere else. 
The reader should keep in mind that these cluster data do not constrain the IFMR at masses 
lower than about 1.7 Mq. As a consequence, the present analytic relations, particularly the 
piecewise linear fit, should not be extrapolated towards lower initial masses. 

We could attempt more complicated fits to the data, but in our opinion the uncertainties 
in mass determinations, and possibly the existence of an intrinsic spread in the IFMR to be 
discussed below, render such efforts probably unnecessary. In fact, we have found polynomial 
fits beyond second order do not produce any relevant improvement in the goodness of fit. 

As an additio nal check, we p e rform ed a direct comparison of our fits with the polynomial 



relationship from iFerrario et al.l (120051 ). that was derived by combining an assumed star 
formation rate and initial mass function for the Galactic Disk, with the observed field WD 
mass distribution. For masses between 2.5 and 6.5 Mq - the range of initial masses covered 
in their study - at a given Mi, differences in the M{ values are less than 0.04 Mq, well below 
the dispersions given above. 

Another point that we wish t o address is the presence of a possible intrinsic spread in 



the reference IFMR, as noticed by iReicu (119961 ) in an analysis of Praesepe WD population. 



With intrinsic spread we mean a signature of diff erential mass loss amo ng stars with the 



same mass and same initial chemical composition. ICatalan et al.l (J2008al ) have investigated 
very recently the spread around their mean semi-empirical IFMR (determined from cluster 
WDs and WDs in common proper motion pairs) and found no clear cut indication that is 
either due to a spread in mass loss or to some other individual stellar properties (e.g. the 
thickness of the envelope layers) that vary between objects. 

To investigate this issue with our data, we made the following test. We considered 
objects in a single cluster, namely Praesepe, with M; values in agreement within their re- 
spective la error bars. There are 8 stars with Mj between 3.045 and 3.543 Mq satisfying 
this constraint, and we assume that - within the present error bars - they share the same 
Mi, equal to the average value (Mi) = 3.207 Mq. For these objects, the average WD mass 
is (Mf) = 0.769 Mq with la rms equal to O.O45M0. To test whether this observed disper- 
sion in Mf is due exclusively to the errors in the Mf estimates, we built a synthetic sample 
by drawing randomly 10000 Mf values for each of the 8 objects, according to a Gaussian 
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distribution centred on (Mf) = 0.769 Mq, with a la spread given by their individual error 
bars. The synthetic distribution has obviously (Mf) = 0.769 Mq, but a Icr rms equal only 
to 0.028Mq. We compared statistically the a values of the two distributions by means of an 
F-test, that returned a probability of P=96.3% that a of the synthetic sample is different 
from (smaller than) the observed one. By fixing, as customary, a threshold of P=95% to 
accept the existence of a real difference between the two a values, our simple test seems 
to point to the existence of a spread in the IFMR, at least in this cluster. Adding 5 more 
objects in the Hyades and in NGC3532 (the [Fe/H] estimates for these two clusters plus 
Praesepe span a range of only ~0.10 dex, and overlap within their associated ~ 2cr error 
bars, see Table [3]) with initial masses within the same overlapping Mi range and repeating 
the previous analysis does not change this conclusion. In fact, for this enlarged sample, 
(Mf) = 0.729 Mq, with a la rms equal to O.IOMq. The probability that the corresponding 
synthetic sample - constructed in the same way as for the 8 Hyades WDs - has a different 
a than observed is P=96.9%. 

However, this evidence for an intrinsic spread (i.e. due to a significantly varying mass 
loss at a fixed mass and chemical composition) in the IFMR is weakened if we take into 
account theoretical expectations for near-solar metallicities, displayed in Fig. [3J Although 
the Mi values estimated for these 13 WDs (or just the 8 Praesepe objects discussed above) 
overlap within the la error bars, the same error bars suggest that they may also cover a 
range of values. Just considering the range covered by the central Mi estimates for these 
13 WDs, theoretical IFMRs (see also another completely inde pendent theoretical IFMR for 



solar metallicity in Fig. 3 of iPrada Moroni fc Stranierd 120071 ) generally display a steep rise 



of Mf with Mi. Over the Mi range between 3.05 and 3.5 Mq, the theoretical IFMRs of 
Fig. [3] predict an increase in Mf by about 0.1 Mq, which corresponds to ~ 2a of the observed 
Mf spread in the 8 Praesepe objects, and to ~ la of the observed spread in the 13-object 
enlarged sample. Without accounting for errors in the estimate of Mf , the predicted change 
within the Mi range allowed by the error bars of the estimated initial masses can potentially 
explain the observed Mf spread. This would imply that the observed spread is not the result 
of star-to-star mass loss variations, but due to, as predicted from theory, the steep increase 
of Mf with Mi in the relevant mass range. 

Moving this analysis to different Mi ranges can hardly provide additional constraints 
on the existence of an intrinsic spread in the IFMR, because of either smaller numbers of 
objects in comparable Mi intervals (when moving to lower Mi values) or a combination of 
much larger error bars on Mi and small numbers (when moving to higher Mi values). It 
may however be interesting to notice the situation at Mi ~ 1.8 Mq, where we have 5 objects 
(belonging to NGC7789 and NGC6819) with formally the same Mi, small errors in both Mi 
and Mf and theoretical predictions for an IFMR essentially flat (see Fig. [3]). By applying the 
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same statistical analysis described before we do not find any statistically significant spread 
in the IFMR for this sample. 

Before concluding this section we discuss briefiy another test of our reference IFMR. In 
the case of NGC2099 we have co nsidered the alterna tive value [Fe /H]=— 0.20 and reduce d 



reddening (E(B-V)=0.23) used by lKalirai et al.l (120051 ) and given by lDeliyannis et al.l (120021 ). 



Note however that Deliyannis (private communication) considers this as a possible but as yet 
unconfirmed value for [Fe/H]. Using this alternative metallicity and reddening we derive an 
age for NGC2099 of 550 Myr which, compared to 320 Myr derived with [Fe/H]=0.09, has a 
dramatic effect in the inferred initial masses (Fig. [7j). Results for NGC2099 stars (#27-#38) 
are repeated in the last part of Table [6] but in this case corresponding to the low [Fe/H] . As 
expected from the relatively large change in the inferred cluster age, changes in Mj due to 
the different cluster parameters generally exceed the uncertainties in Mi determination. 

This test, based on estimates that are admittedly not confirmed, is just to exemplify in 
a quantitative way how crucial the assessment of the cluster [Fe/H] and E(B-V) values is. 

We will return to this issue in Sect. 4. 



3. Global error-bar estimates 

The method used to determine the IFMR has been described in Section [2], where we 
also presented our reference IFMR. We now focus on a detailed determination of the global 
uncertainties both in the initial and final masses for each star in our sample. Since the global 
errors have contributions from various sources, they have been estimated by means of Monte 
Carlo simulations. In the following we detail the procedure followed for each star in our 
sample. 



3.1. White Dwarf Uncertainties 

We have tried to cover all major sources of uncertainties in the WD mass and cooling 
age determinations, that are going to be discussed in this section. While we have used 
SOO models as our standard set of WD models, all the WD evolutionary sequences to test 
the effect of WD input physics and chemical stratification have been computed with the 
LPCODE described below. This should not invalidate our procedure, since we are using 
LPCODE models only to determine differential changes in the WD properties (e.g. core 
composition), not absolute values. In what follows, X represents the logarithmic values of 
either the WD mass or the cooling age. 
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- Observational uncertainties. The central values and la uncertainties for log*? and Tcs 
are listed in Table [H We leave them unmodified as they are given in the original papers 
quoted in Sect. 2.1. 

- Systematic uncertainties from WD cooling tracks. In addition to our standard choice 
of WD cooling tracks (SOO), we have specifically computed a completely independent set of 
WD cooling models with the LPCODE. The aim is to estimate the systematic uncertainties 
in the inferred WD properties that arises from using different WD models. We define the 
la uncertainty by 

0"syst(-^) = I^SOO ~ -^LP|; (3) 

where Xsoo is derived from SOO models and Xlp from models c omputed with the LPCODE . 
WD models computed with the LPCODE have been used in ISerenelli fc Fukugital (120071 ). 
where the evolution of a set of stellar models from 1 to 8 Mq and solar metallicity has 
been consistently followed from the main sequence phase up to the WD phase, including 
the thermal pulse (TP) AGB phase. It should be noted that, in addition to differences in 
numerics, CTgystl-'^) also accounts for differences in input physics between the two sets of 
m odels. While S OO and LPCODE models have used th e same neutrino energy loss rates 



(lltoh et al.lll996l ). conductive opacities in the CO cores (jltoh et al.lll993l ) and very similar 
hydrogen envelope thickness and core composition, they differ i n other relevant asp ects 
such as the equat ion of state and low temperature opacities (see lAlthaus et al.l l2003l and 
Salaris et al.l l2000l for details about the LPCODE and SOO input physics respectively) and 
the inclusion of gravitational settling in the model envelopes (SOO models mimic the effects 
of gravitational settling by adopting a pure hydrogen envelope and fixed chemical H/He 
and He/C/0 interfaces while in LPCODE the chemical profiles are evolved according to 
time-dependent equations of element diffusion). 

- Neutrino energy loss rates. For the WD evolutionary sequences we are considering in 
this work, plasma neutrinos always contribute at least 90% of the total neutrino emission. 
Uncertainties in the theoretical calculation of plasma neutrino emission rates, within the 
framework of the standard model of particle physics, are at most a few percent (Fukugita, 
private communication) and thus we do not expect them to represent an important source 
of uncertainty in determining WD cooling ages. Several additional cooling mechanisms 
that could operate in WD interiors have been proposed, e.g. enhanced neutrino cooling 
due to a non-null neutrino magnetic dipole m oment, axion production or other particles 
like Weakly Interacting Massive Particles (see iRaffeltl Il996l for a complete discussion of 



these mechanisms). It is not our aim here to consider the effects on the WD properties - 
particularly on the cooling rate - of all such mechanisms, especially because they remain 
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purely hypotheticao We do not want, however, ignore the possibihty that some additional 
coohng might be present. Consequently, we have adopted a more conservative point of view 
and have computed additional W D cooling models where standard neutrino energy loss rates 
(JHaft et al.lll994l : lltoh et al.lll996l ) have been multiplied by factors of 2 and 0.5 respectively to 
take into account the possibility of such additional cooling mechanisms. These large factors 
were chosen to represent our conservative approach. We define the WD mass (and cooling 
age) Icr uncertainties due to neutrino energy losses as 



a.(X) 



\Xo — Xi 



0.51 



(4) 



where the subindices refer to quantities corresponding to the WD cooling models with the 
modified neutrino emission rates as described above. 

- Conductive opacity. For the physical conditions and composition of a representative 
set of WD m odels, we have c omputed conductive opacities according to both Itoh and 
collaborators (jltoh et al.l (119931 ) and references therein - our refe rences set of conductive 
opacities) and Potekhin and collaborators (see ICassisi et al.l 120071 for the most up to date 
calculations). We find differences between both sets of calculations to be usually below 20% 
for the range of temperatures and densities relevant to the WD models used in this paper. In 
analogy with the neutrino emission rates, we have estimated the effect of electron conduction 
uncertainties by computing additional WD models where the conductive opacities from Itoh 
and collaborators have been multiplied by factors of 1.25 and 0.8 and have used these models 
to define the WD mass and cooling age la uncertainties due to errors in conductive opacities 
calculations as 

1-^1.25 ^ -^0.8 1 



CT.iX) 



(5) 



- WD core composition. The mass fraction of the central oxygen Xq in SOO WD models 
ranges from 0.65 up to 0.80 depending on the mass of the WD. The internal CO stratification 
is however subject to uncertainties due to the uncertainty on the C^^ + a reaction rate (see, 
e.g. SOO) and, very important, the treatment of mixing (semiconvection, b reathing pulses 



and the mechanism to suppress them) during the central He-burning phase (IStraniero et al. 



^ Among: the non-stan dard cooling me c hanism s, axions are probably the best motivated. Recently, 
BischofF-Kim et al.l (|2008l ) and llsern et al.l (J2008f ) have used WD pulsation pro perties and the lum i nosit: 



function of local Disk WDs, respe ctively, to constrai n the axion mass iria- While iBischoff-Kim et al 



(200. 



find an upper limit nia = 13 meV, llsern et al.l ( 20081 ) rule out values as low as 10 meV and give a preferred 
value of nia ~ 5 meV. With this mass value, axion cooling is never the dominating WD cooling mecha- 
nism (contrary to neutrino cooling, that dominates for hot WDs) but it is strong enough that it can have a 
moderate effect on WD age determinations for WDs with luminosities in the range —1.2 < logi/LQ < —1.7. 
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20031 ). These effects are larger than the expected change in CO profiles due to the small 
metallicity range covered by the cluster sample. In our analysis we have been conservative 
and have defined the la uncertainty due uncertainties in the core composition as 

acore(X) = , (6) 

where Xo.g and X0.3 are derived from WD cooling tracks where the original chemical profiles 
have been rescaled to have Xo= 0.9 and 0.3 at the center, respectively. 



- Hydrogen envelope thickness (Mh). As shown by, e.g.. [Provencal et al.l (Il998l ) the H- 
envelope thickness of DA WDs displays a range of values. To include this effect in our error 
analysis, in addition to the standard thick envelopes considered in our LPCODE calculations 
(usually 10~^ M©, but larger values were used for WD masses below 0.6 M©) we have 
computed also cooling sequences with thinner H-layers of 10~^ Mq. Similar to the above 
procedures, we have defined the corresponding la errors 

\ Y Y I 

ac„v(A) = ; (7) 

where Xthick and Xthin are derived from WD models with thick and thin hydrogen envelopes 
respec tively. We have perfo rmed an additional test to explore very thin hydrogen envelopes 



as m 



Catalan et al.l l2008al ) by computing a 0.6 and a 1 M© WD models with hydrogen 
envelopes of Mh ~ 10"^'^ M0. When compared to models with thick envelopes, we have 
found that differences in the cooling ages never exceeded 11% in both cases and are usually 
at the level of 7% for the 0.60 Mq model and 5% for the more massive one in the range 
of effective temperatures of interest to this work. Asteroseismology studies of ZZ Ceti stars 
disclose a wide rang e of hydrogen envelope thi ckness, ranging between ~ 10~^ M© and 



10 ^ Mq (see, e.g.. lCastanheira fc Keplerll2008l ). In this regard, the changes in the cooling 



ages due to uncertainties in the H-envelope thickness quoted above can probably be taken 
as robust upper limits. 

For neutrino emission rates, conductive opacities, core composition and hydrogen-envelope 
thickness we have assumed that uncertainties distribute normally. In the case of systematic 
uncertainties from different WD cooling tracks we have assumed uniform distributions of 
uncertainties (/x) which are defined, in each case, by J_^^^^^ fxdX = 0.683. We tabulate 
a,y{X), cr^{X), acorc(Ar), cr^rwiX) and asyst(^) in (log5f,Teff) grids for using them in the Monte 
Carlo simulations. We summarize the results for the WD uncertainties, both for mass and 
cooling age, in Table HJ where for a small subset of (log (7, Teg) values, we give the 1 — a 
fractional uncertainties in WD mass and cooling age as defined by Eqs. [3113 but transformed 
to linear scale to facilitate interpretation of the results. 
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3.2. Progenitor Star Uncertainties 

Progenitor masses Mi are obtained from the estimated progenitor lifetime and stellar 
models. In addition to the uncertainty sources discussed in the previous section for tcooi, the 
determination of Mj is affected by: 

- Cluster age. Progenitor lifetimes are obtained simply as tprog = tdns — ^cooi, so un- 
certainties in tcius directly affect the determination of Mi. Central values and adopted la 
uncertainties for the cluster ages determined for this paper are listed in Table [31 As our 
standard choice, we have adopted uniform distributions ftitdus) for cluster age uncertainties, 
where ft is a constant defined by the simple relation /_J"''"° ftdtdus = 0.683. 

- Stellar metallicity (A[Fe/H]). A[Fe/H] affects the determination of Mj in two dif- 
ferent ways. The first and most important relates to the determination of the cluster age 
(mainly through the change in the distance modulus obtained via MS-fitting) and has al- 
ready been taken into account in the uncertainty of tcius- Additionally, A[Fe/H] introduces 
uncertainty in the determination of the progenitor mass because, for different metallicities, 
a fixed tprog corresponds to different progenitor masses. We treat this effect by adopting 
the central values and uncertainties in the metallicities listed in Table [3] to generate stellar 
metallicity distributions (normal distributions are our standard choice). For a fixed tprog we 
then compute M^ for the corresponding metallicity distribution. 

- Different stellar models and isochrones. In order to estimate the systematic uncer- 
tainties arising from using different sets of stellar models and isochrones we have proceeded 
in a similar fashion as in the case of the WD cooling tracks and define as the la uncertainty 

Crsyst(Mi) = |Mi3aSTI " Mi^Paduah (8) 

where Mi BaSTi represents the progenitor masses derived with Basti stellar models and isochrones 
and Mi Padua those derived with Padua models and isochrones. We have adopted a uniform 
distribution for systematic uncertainties. It should be noted here again that consistency in 
the calculations requires the same set of models and isochrones be used to derive Mi and 
tprog, i.e. the cluster age. It is in principle not correct to adopt a cluster age derived with, 
let us say the BaSTI isochrones, and then the Padua stellar models to derive Mi from the 
calculated tprog. As already stated in the Introduction, such inconsistencies are ubiquitous 
in the current literature and can introduce noticeable systematic changes in the derived Mi, 
particularly for Mi ^ SM©. 

We close this section by stressing that both codes, BaSTI and Padua, are employing up- 
to-date physics input, and therefore are, at least nominally, quite similar in their treatment of 
stellar evolution. Systematic uncertainties arising from code-to-code differences are therefore 
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possibly underestimated in the above error analysis, if they arise from different physical 
assumptions, as, for example, the inclusion or neglect of convective overshooting, which we 
did not include in the discussion of this section (we provide separate results for models 
without core overshooting). 

When considering models without core overshooting, the systematic effect of using dif- 
ferent models and isochrones was taken into account with a slightly different approach be- 
cause Padua calculations are only available for one chemical composition, i.e. Z=0.019 
and y=0.273. Th i s pair of {Y, Z) values corresponds to [Fe/H]=0.04 when considering the 



Grevesse fc Noelsl ( 119931 ) solar mixture (that provides Z/X=0.0245±0.005 for the present 



Sun) e mployed in those stellar model calculations. Also BaSTI models employ the lGrevesse fc Noels 



( 119931 ) solar mixture, and the {Y, Z) pair closest to the Padua one around the solar metal- 
licity is Z=0.0198, Y =0.2734, that corresponds to [Fe/H]=0.06. In this case, and just for 
estimating this systematic uncertainty, we have computed all progenitor masses for both 
BaSTI and Padua non-overshooting models assuming a fixed [Fe/H]=0.04 for Padua and 
[Fe/H]=0.06 for BaSTI. As above, the difference between the masses has been adopted as a 
measure of the systematic uncertainty and we have assumed this difference does not depend 
on the actual [Fe/H] value. This is in fact a reasonable approximation and, in addition, 
systematic uncertainties from isochrones and models turn out to be a small contribution to 
the overall uncertainty budget, so our simplified treatment should not affect our results and 
conclusions. 



3.3. Monte Carlo simulations 

For each star in our sample, the Monte Carlo (MC) simulations consist of the following 
steps: 

- generate {logg-^, Teffjjjl^ distributions according to observational uncertainties and use 
SOO WD models to get the initial WD mass {M°j}|l^ and cooling age {t°ooi jljl^ distributions. 
Here A'^(= 10^) is the total number of trials in the MC simulations and quantities with the 
j suffix denote results of each one of the trials; 

- For each trial j include additional WD mass and cooling age uncertainties according 
to the details given in Section [XT! Each mass in the the final WD mass distribution {Mf j}|l^ 
(analogously for the WD cooling age) is given by 

Mf,j = Mf. + J2 ^^J,k, (9) 

k=l,5 

where each (5Mj k, the individual contributions to the WD mass global uncertainty, is ran- 
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domly drawn from distributions constructed as described in Section I3.1[ The final WD 
cooling age distribution {tcooijjjli i^ computed in the same way; 

- generate the cluster age distribution {tcius,j}jli ^^'^ then the progenitor age distribution 
{if^progjljli, where for each trial tprogj = t^usj - ^cooij; 

- From {tprogjjjli get the initial progenitor mass distribution {Mjjjjl^ using the BaSTI 
stellar models; 

- The final progenitor mass distribution {Mijjjl^ is obtained by adding the star metal- 
licity and stellar models systematic uncertainties, i.e. Mij = M°j + (5Mi [pe/H] + 5Mi_syst as 
discussed in Section I 



4. Results and comparisons with theory 

4.1. Global error assessment 

Table displays the results about cooling age and mass for all the WDs in our sample. 
The central values have been obtained using SOO WD models and the uncertainties are the 
result of the MC simulations described in the previous section and thus include all sources 
of uncertainty considered in this work. We also provide, as additional information, the 
bolometric luminosity of the WDs, as obtained from the interpolation among the WD model 
grid in log f? and Tgfj. It is important to notice that none of the objects have yet entered the 
crystallization phase. This minimizes the effect of uncertainties in the interior CO profiles, 
that play a very important role during the phase separation associated to the crystallization 
process (SOO). 

Our reference IFMR relation (BaSTI with core overshooting plus SOO models) is dis- 
played in Table [6] (labelled "ov"). For the WD masses, the quoted uncertainties are the 
average of cr_(Mf) and a+{M{) given in Table [51 For the progenitor masses we give the re- 
sults of our Monte Carlo simulations, that show strongly asymmetric error bars arising from 
the highly non-linear relation between evolutionary lifetime and mass for MS stars. Results 
for BaSTI models without core overshooting are also shown (labelled "non-ov"). 

As already discussed, our estimates of the global uncertainties in the IFMR include 
systematic effects due to the use of different sets of stellar models and isochrones, and 
different WD cooling models. It is instructive, however, to try to single out the influence of 
the stellar models on the final error. To this aim. Fig. H] shows a comparison between the 
progenitor masses obtained with two different sets of stellar models and isochrones (BaSTI 
and Padua). In both cases we used SOO WD models. When compared with the error bars 
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on Mi shown in Fig. [3l Fig. H] leads us to conclude that systematic uncertainties from using 
different sets of isochrones and stellar models do not contribute appreciably to the global 
error budget on Mj, as long as the assumptions concerning the stellar physics are similar. We 
emphasize again that to achieve internal consistency the same set of isochrones and stellar 
models must be used to determine the cluster ages and to obtain progenitor masses from the 
progenitor lifetimes. 

We next study the changes in the IFMR from using different WD models. Figure 
displays our reference IFMR (SOO + BaSTI with overshooting) and that resulting from using 
the LPCODE WD models. Although WD masses derived from both sets of WD models are 
very similar, the effect on the progenitor mass is much larger, especially for progenitor masses 
above ~5 Mq. The changes are due to differences in the WD cooling ages obtained with the 
different WD models. Both SOO and LPCODE models give WD ages consistent within the 
observational errors in \ogg and T^e-, however the short lifetime of massive stars makes the 
determination of Mi very sensitive to small differences in tcooi for higher progenitor masses. 

Figure E] gives an overview of the relative contributions of observational and model 
(cooling tracks, including systematic effects due to variations of input physics and chemical 
stratification) uncertainties to the total error budget for the final WD masses and ages. We 
display also the relative contribution of cluster age (due to isochrones, cluster metallicity, 
reddening, distance) and WD age (due to cooling tracks) uncertainties to the total error 
budget for the progenitor masses and ages. These fractional contributions are expressed 
in terms of (o"i/crtot)^ for each WD, denoted by its number (as in Tabled]), whilst vertical 
dashed lines separate the various clusters. The uncertainty of the WD mass appears to be 
completely dominated by the observational errors in log g and Teg , while in case of WD ages 
the systematic uncertainties of the cooling tracks dominate for some clusters (NGC3532 and 
NGC2099). For the progenitor ages and masses (lower panels of Fig. E]) cluster ages are the 
largest contributors to the total error, except mainly for NGC2099. 

In closing this section, we comment on the possibility that systematic uncertainties (to 
be added to the q uoted errors) may be p resent in the empirical determination of WDs Teg 



and log5f(see, e.g., iNapiwotzki et al.lll999l ). We have not included this source of uncertainty 
in our global error assessment simply because it is not quantified in the literature. We 
have performed, however, a simple test to understand how such hypothetical uncertainties 
may influence the determinations of WD masses and cooling ages. We have considered 
that Teff determinations have systematic uncertainties amounting up to ±1000 K, that we 
represent by means of a uniform distribution. Analogously, in the case of \ogg we have 
assumed an additional systematic uncertainty of ±0.15 dex. Note that these are probably 
very generous systematic errors, because we assume all WDs are similarly affected. Our 
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results show that uncertainties in logg have a more noticeable effect than those on Tes- This 
is because WD masses are more robustly predicted by WD models than cooling ages, and 
consequently the error budget is dominated by observational uncertainties (as already shown 
in Fig. E]). For about 40% of our sample, uncertainties in WD masses are about a factor of 
2 larger when systematic uncertainties as described above are included; these are the stars 
that show the smallest observational uncertainties in Table [H For the rest of the stars in 
our sample, uncertainties are increased by not more than a factor of 1.5. Cooling ages, on 
the other hand, are much less affected by systematic uncertainties, because the total error 
budget here is heavily influenced by input physics and systematics in WD models. As a 
consequence, initial mass determinations are only very mildly affected. Indeed, the largest 
increase in positive uncertainties for Mi is only 50%, with only 4 objects having changes 
above 40%. There are 40 objects for which positive uncertainties increase by less than 
25%. Negative uncertainties, which are always smaller, are less affected, and for 49 objects, 
increases are smaller than 25%. We thus conclude that adding systematic uncertainties 
to WD observational parameters would only have a very modest influence on the inferred 
IFMR. Based on these results and also considering that any attempt to quantify systematic 
uncertainties remains, at this point, arbitrary, we have chosen not to include these effects 
neither in our reference IFMR nor in the discussion of our results. 



4.2. Comparison between predicted and derived IFMRs 

Figure [3] compares the inferred IFMR with that predicted by the same BaSTI models 
(we chose to display the BaSTI results for [Fe/H]=0.06) used in the determination of tp^og and 
Mj. We display results obtained from progenitor models and isochrones with and without 
MS core overshooting. 

Cluster ages determined from models without overshooting are smaller (see Table [3]), 
hence tprog is decreased, resulting in larger progenitor masses. 

As a consequence, although the associated error bars are large, models without over- 
shooting predict for some stars progenitors masses well above the maximum possible val- 
ues allowed by theory. Examples are the WD in the Pleiades and two of the WDs in 
NGC2516 (see Table [6]), for which the estimated progenitor masses are between ~ 12 and 
~ 16 Mq. Theoretical stellar evolution models predict an upper limit to the progenitors of 
carbon/oxygen white dwarfs in the range between ~ 6 and ~ 8 Mq, the exact values depend- 
ing on the chemical composition of the models and on the treatment of mixing in the cores; 
lower upper limits are obtained when core overshooting during the central H-burning phase 
is included. In several other cases in Table M even negative progenitor ages were obtained. 
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The relationship between Mi and the core mass at the onset of the thermal pulses 
(MciTp) is also displayed in Fig. [31 For models without core overshooting, the theoretical 
IFMR lies largely above the semi-empirical determination, especially for M^ < SMq, while 
the theoretical M; — MdTp relationship is in better agreement with the data. This would 
imply that the progenitor stars did not experience any thermal pulses at all. 

On the other hand, for core overshooting models the theoretical IFMR follows much 
better the semi-empirical results, and no very high initial masses are typically derived in this 
case. The theoretical M[ — Mcitp relation constitutes a lower envelope to the data, consistent 
with general expectations about AGB evolution. This would lead us to conclude that core 
overshooting during the MS results in a better agreement with the derived IFMR, because 
of the internal consistency about the end-product of stars with inferred Mj. 

Before moving forward in our comparison with theory, we wish to comment briefly on 
NGC2099, that represents an interesting case and a warning about the role played by the 
cluster [Fe/H] and E(B-V) estimates. As already mentioned, we have considered two possible 
values for its metallicity and reddening. Results for both choices are given in Tableland are 
also plotted in Fig. [71 It becomes apparent from the figure that data points fall systematically 
below even the theoretical Mi — MdTP (models and IFMR with overshooting) relation for 
[Fe/H]= 0.09. The same result is obtained if models without core overshooting are used 
ins tead. This in c onsist ency is strongly alleviated if the alternative [Fe/H] and E(B-V) used 



by iKalirai et al.l (l2005l . see Table [3]) are adopted. While the aim of this paper is not to use 
the IFMR to improve cluster properties, this example demonstrates how informative it can 
be, when systematic effects are well controlled. 

Focusing again on the core overshooting IFMR (Fig. [3l top panel) we can finally try 
to obtain some constraints on TP-AGB evolution by comparing the inferred IFMR with 
that resulting from sequences with different TP-AGB modelling. To do this, we first com- 
pare in Fig. [3l the Mi — MdTP relations from BaSTI and LPCODE stellar models (thick 
dotted and thin dashed lines respectively). Both relations are very similar, regardless of 
differences in the treatment of convective and additional mixing episodes. We recall that in 
BaSTI instantaneous mixing in the overshooting region beyond the Schwarzschild boundary 
of the H-burning convective core is ass umed, as well as semiconv ective mixing during He-core 



burning, which is briefly described in [Pietrinferni et al.l (120041 ). where also the appropriate 



references are given. On the other hand, no extra mixing is taken into account at the bottom 
of convective envelopes. LPCODE models include diffusive mixing due to overshooting at 
all convective boundaries. We also point out that evolutionary timescales until the onset 
of the TPs are very similar in BaSTI models with overshooting and the LPCODE models. 
Since TP-AGB lifetimes are very short and their contribution to the the total progenitor 
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age can be neglected, we can safely use semi-empirical data obtained from BaSTI models for 
comparisons with the LPCODE theoretical IFMR. 

Having assessed the similarity of the M^ — Mcyvp relations, differences in the IFMRs 
predicted by BaSTI and LPCODE models can be traced to differences during the TP- 
AGB phase. While BaSTI models do not include extra mixing during the thermal pulses, 
LPCODE includes, as mentioned above, diffusive overshooting at all convective boundaries. 
The theoretical IFMRs are also shown in Fig.Owith thick (BaSTI) and thin (LPCODE) solid 
lines and do indeed differ. This is unrelated to the choice of mass loss during the AGB phase, 
and it is mainly due to the presence of very strong third dredge-up episodes in LPCODE 
sequences, which inhibit the growth of the H-free core for models with Mi ^ SM©. These 
strong third dredge-up episodes can be traced back mainly to the presence of overs hooting 



at the base of the thermal pulse driven convective zone. This was demonstrated by iHerwig 



(J2OOOI ) and confirmed by us running test models where convective overshooting was included 
only at separately selected convective boundaries during the AGB evolution. Due to strong 
third dredge-up events, the final mass of LPCODE sequences with Mi ^ 3Mq is very similar 
to MciTP, the core mass at the first thermal pulse (Fig. [3], top panel). On the other hand 
in BaSTI sequences the final mass is mainly set by the length of the TP-AGB phase and, 
thus, by the adopted mass loss rate. The apparent failure of the theoretical LPCODE 
IFMR to match the final masses in the range 2.75Mq <^ M^ ^ S.SM© is an indication that 
such strong dredge-up episodes do not take place in real stars, and therefore the adopted 
exponentially decaying overshooting at the base of the thermal pulse driven convective zone 
in the LPCODE models is too stro ng. This allows us to conclude that the overshooting 



parameter / (see iHerwig et al.lll997l . for a definition of /) must be significantly smaller than 
/ = 0.016 employed at the convective core of upper MS stars. This is in line with recent 
hydrodynamical simulations of AGB thermal pul ses which point to a value of / ^ 0.01 at 



the bottom of the pulse driven convective zone (IHerwig et al.l 120071 ). On the other hand, 
given the large error bars of the data points, the existence of strong third dredge-up episodes 
for stars with Mi ^ 5Mq cannot be excluded by the same arguments. The IFMR may 
therefore teach us that overshooting varies a lot between different convective layers and may 
also depend on stellar mass. 



5. Summary and conclusions 

Since its earliest determination, the semi-empirical IFMR has indicated a monotonic 
rise of the final WD mass with the initial mass of the progenitor. It is also strong evi- 
dence that single stars of intermediate mass do not develop degenerate cores close to the 
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Chandrasekhar limit, most likely due to the effect of strong mass loss on the AGB. These 
qualitative results have remained unchanged in spite of the significant improvement both in 
the amount and accuracy of observational data and in the theoretical models. Nevertheless 
there are still many open questions, such as a dependence on metallicity, the slope of the 
IFMR at the upper (initial) mass end, changes in slope, and spread around a mean relation. 
To answer these questions, not only increasingly accurate observations are needed, but also 
the systematic uncertainties coming from the theoretical input - which is necessary to es- 
tablish the semi-empirical IFMR - have to be controlled. These are related to evolutionary 
tracks and isochrones that determine the total age of the observed WDs and cooling tracks 
used for deriving the WD age. Our paper was intended to investigate the influence of such 
systematic uncertainties and to derive robust features of the IFMR. We also investigated 
what uncertainties in the theory of stellar evolution would lead to inconsistent results, and 
therefore could be restricted or excluded. 

From the available WD data we selected a total of 52 objects in 10 open clusters, plus 
Sirius B. There are, in fact, even more objects available, but we disregarded those for which 
our cluster distance determination (Sect. 2) would not be possible because of the lack of 
calibrating field stars at that r netallicity (NGC6791 a t [Fe/H] = 0.45). We also did not take 
into account the objects from ICatalan et al.l (l2008bl ). because of the more indirect way to 



obtain the age of the mostly unevolved binary companion. Inclusion of such objects would 
weaken our aim of highest possible self-consistency. 

For the deter mination of mass an d cooling age of th e WDs we used diff erent sets of 



cooling track, from ISalaris et al.l (120001 ) and the LPCODE (jAlthaus et al.l 120031 ) to evaluate 
the influence of different codes, input physics and chemical stratiflcation (Sect. 3.1). We 
estimated the error ranges quite conservatively and to evaluate the global error budget we 
used Monte Carlo simulations (Sect. 3.3), for which we assumed either a normal or uniform 
distribution of the different error sources. We included the effect of observational errors, too. 

Results concerning cooling ages and final masses are quite robust: the accuracy of 
the WD masses depend completely on the observational uncertainties, while their ages is 
infiuenced in most cases (in terms of number of clusters) by the uncertainties of the cooling 
tracks, mostly systematic differences between different sets of cooling tracks and changes in 
the neutrino emission rates. 

This result implies that any spread in Mf , as seen in Fig. [21 which is outside the error 
bar is real, in the sense that Mf does vary at the same (within the errors in M;) initial mass. 
Whether this spread is consistent with theoretical expectations, related to a steep increase 
of Mf with Mi in that mass range, or whether it is due to star-to-star variations in the total 
mass lost, cannot be decided yet, as the error bars in Mj are too large even in the case of 
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the Praesepe WDs, which we discussed in detail in Sect. 2.4. 

Figure \5\ demonstrates that uncertainties due to different WD codes with similar, but 
not identical physics and model details are essentially irrelevant for the estimate of Mf. 
However, differences in the cooling ages have a significant influence on Mi, since they af- 
fect directly progenitor lifetimes, and the more massive progenitors evolve on very short 
timescales. Therefore the two different sets of cooling tracks lead to quite different values for 
Mi, particularly for initial masses above 5 Mq. In this regime, M[ differences are between ~ 
0.5 and 2.0 Mq. 

The determination of Mi is more complex, as it involves several steps. The age deter- 
mination of the cluster or binary system requires metallicity, reddening and distance. We 
have used literature values for the former two, and our own homogenous determination for 
the latter. Uncertainties from both steps have been taken into account. Using appropriate 
isochrones the age is then obtained, which, together with the WD cooling age, results in the 
pre- WD lifetime. As this is, to an accuracy of 1% or slightly more, identical to the lifetime 
from the MS to the end of the central He-burning phase, all uncertain ties concerning the AGB 



evolution are irrelev ant. We employed different sets of stellar models (jPietrinferni et al.ll2004 



Girardi et al.ll2000l ). which differ slightly in the detailed treatment of the input physics (but 
both include core overshooting during the MS) as well as additional models from the former 
source (BaSTI), which do not consider convective overshooting. The first effect - different 
tracks - has a very minor effect on Mi, as we showed in Fig. HI Note that this includes both 
the cluster age determination and the relation between progenitor age and mass. The neglect 
of overshooting, however, leads to internal inconsistencies in the semi-empirical IFMR: at 
the high end, initial masses are predicted, for which the models themselves predict that no 
electron degenerate CO- WD would result from the evolution. By discriminating between the 
different error sources (but omitting the overshooting effect). Fig. demonstrates that the 
cluster age uncertainties (due to uncertainties in the cluster [Fe/H] and reddening estimates) 
dominate the error on Mi, although the above-mentioned WD age can be significant for the 
more massive progenitors, too. 

The case of NGC2099 illustrates the importance of accurate metallicities and reddenings: 
Figure [7] shows how drastically the semi-empirical IFMR changes, two different estimates of 
[Fe/H] and E{B — V) are assumed. If one accepts this uncertainty as being real, NGC2099 
cannot be used at all to learn about the IFMR. 

After we had determined the full extent of the uncertainty associated to our semi- 
empirical IFMR, we compared it with the predictions from the theoretical calculations, 
both from BaSTI and LPCODE. Both differ in their treatment of convective overshooting 
during the AGB phase. One should recall that the BaSTI models treat that phase with 
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synthetic models, while LPCODE follows the evolution in full detail. Figure [3] shows them 
in comparison with the semi-empirical IFMR and the relation between M^ and the mass 
of the CO-core at the first thermal pulse, Mcitp- As long as stars do not avoid the AGB 
phase, Mf > M^tp should be fulfilled. Therefore, this relation is a strong lower limit 
to the semi-empirical IFMR. The agreement between the Mi — MdTP relations from both 
codes is excellent; there are no significant systematic uncertainties in its prediction. This 
again emphasizes the need for convective overshooting during the MS phase: without it a 
significantly larger numbers of objects falls below this limiting line. 

On the other hand, the theoretical IFMR from the LPCODE tracks is not consistent with 
the data: In these models overshooting during the AGB, mainly from the lower boundary of 
the pulse-driven convective zone, limits the growth of the CO-core such that the final core size 
is almost identical to that at the first pulse. However, most objects have final masses clearly 
above this line, such that we concluded that the extent of overshooting during the AGB phase 
must be significantly smaller than on the MS. This is in agreement with hydrodynamical 
simulations. 

The theoretical IFMRs displayed in Figure [31 and in particular those resulting from the 
BaSTI models display the largest gradient in the same mass range where a significant spread 
in Mf is observed. We demonstrated that indeed that spread could arise solely from this 
local steep gradient, such that no other explanation would be needed. If it is not, one would 
have to invoke a spread in the total mass lost in a given, small mass range. This is not 
exceptional, as the morphology of horizontal branches in Galactic globular cluster shows, 
but the reasons for such a star-to-star variation are unclear. In other mass ranges the total 
error budget is too large to discuss the reality and reason for apparent spreads in Mf. 

Figure [2] presents our reference semi-empirical IFMR, obtained using BaSTI isochrones 
and tracks with overshooting an d SOO cooling models . The error bars include all errors we 



have investigated. Compared to iFerrario et al.l (120051 ) our error bars are necessarily larger 



(we consider more sources of error), but on average only by less than a factor of two, and 
mostly for Mj. We also present a simple linear fit through all the data as well as a piecewise 
fit, which reflects the properties of the theoretical IFMR based on the BaSTI models. 

In summary, our most important results are 

1. None of the WDs employed in current IFMR determinations is close to the Chan- 
drasekhar mass, not even the progeny of the more massive intermediate mass stars; 

2. stellar models without convective overshooting during core hydrogen burning lead to 
internal inconsistencies in the semi-empirical IFMR; 
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3. overshooting from convective boundaries during the AGB phase must be significantly 
reduced (compared to the case of convective cores along the MS) to reproduce the 
observed Mf values; 

4. the uncertainty in Mf is dominated by observational errors; 

5. the uncertainty in Mi has several reasons: both cluster parameters and isochrone details 
influence the cluster age and thus the progenitor mass; the uncertainty on the WD 
cooling age can sometimes also be the dominant factor; 

6. the observed dispersion in Mf at approximately constant Mj (3.0-3.5 Mq), in particular 
for the Praesepe objects, appears to be real. It may follow from the steep increase of 
Mf with Mi predicted by theoretical IFMRs in that M; range, rather than being caused 
by a spread in mass loss among cluster AGB stars; 

7. from the general agreement between the theoretical IFMRs and the semiempirical 
data, we find no evidence that the mass loss prescriptions used in the stellar evolution 
calculations grossly disagree with the total mass actually lost by low- and intermediate- 
mass stars; 

8. the case of NGC2099 illustrates the necessity for accurate cluster parameters; reliable 
composition data and reddening are required; 

9. as long as up-to-date input physics is used for the stellar models, systematic uncertain- 
ties do not change the overall appearance of the semi-empirical IFMR as determined 
recently by several groups. 

An extension of this analysis to a larger range of [Fe/H] is needed, to investigate a 
possible metallicity dependence of the IFMR, given the small [Fe/H] range spanned by the 
systems analyzed in this study (see Table [3]). In addition to the relevant observational data, 
this requires also a careful study of how to extend to a wider [Fe/H] range the empirical 
MS fitting technique adopted here to determine cluster distances, if we wish to maintain 
the same degree of homogeneity in the derivation of the cluster ages. An extension to lower 
metallicities seems to be preferable, as both calibrating objects and suitable clusters should 
exist, contrary to the case of an extension to higher, that is super-solar clusters. 

As most of the clusters used for deriving the IFMR ha ve nearly solar metall icities, the 



question arises, what the true solar chemical composition is? lAsplund et al.l (120051 ) have pre- 
sented new determinations of the solar heavy element abundance, drastically lower than the 
standard composition assumed here and in the stellar evolution tracks employed. Both BaSTI 
and LPCODE progenitor models consider a standard solar mixture with Z/X ~0.0245, 
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whereas lAsplund et al.l (120051 ) redetermination gives ZjX ~0.0165. Assuming that the dif- 
ferential abundance analyses for our open cluster sample provides the true abundances rel- 
ative to the sun, also the clusters should be less metal-rich than thought. It would be inter- 
esting to see how this will affect the semi-empirical IFMR and its agreement with theoretical 
predictions. An extension of this investigation by using in the computation of progenitor 
models the new solar abundances is currently under way. 
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Table 1. Adopted data for our WD sample, taken from the literature (see Sect. 2.1) 



System 


WD name 


Teff [K] 


cri%s) 


log(^) [cm s 2 


^(log(^)) 


ID 




Pleiades 


LB1497 


32841 


170 


8.630 


0.040 


1 


Hyades 


0352+098 


16630 


350 


8.160 


0.050 


2 




0406+169 


15180 


350 


8.300 


0.050 


3 




0421+162 


19570 


350 


8.090 


0.050 


4 




0425+168 


24420 


350 


8.110 


0.050 


5 




0431+125 


21340 


350 


8.040 


0.050 


6 




0438+108 


27390 


350 


8.070 


0.050 


7 




0437+138 


15340 


350 


8.260 


0.050 


8 


Praesepe 


0836+197 


21950 


350 


8.450 


0.050 


9 




0836+201 


16630 


350 


8.010 


0.050 


10 




0836+199 


14060 


630 


8.340 


0.060 


11 




0837+199 


17100 


350 


8.320 


0.050 


12 




0837+218 


16833 


250 


8.390 


0.030 


13 




0837+185 


14748 


400 


8.240 


0.050 


14 




0840+200 


14180 


350 


8.230 


0.050 


15 




0833+194 


14999 


250 


8.180 


0.040 


16 




0840+190 


14765 


270 


8.210 


0.030 


17 




0840+205 


14527 


390 


8.240 


0.040 


18 




0843+184 


14498 


200 


8.220 


0.040 


19 


NGC2516 


2516-1 


28170 


310 


8.480 


0.170 


20 




2516-2 


34200 


610 


8.600 


0.110 


21 




2516-3 


26870 


330 


8.550 


0.070 


22 




2516-5 


30760 


420 


8.700 


0.120 


23 


NGC3532 


3532-8 


23370 


1065 


7.713 


0.148 


24 




3532-9 


29800 


616 


7.827 


0.229 


25 




3532-10 


19270 


974 


8.143 


0.266 


26 


NGC2099 (M37) 


2099- WD2 


19900 


900 


8.110 


0.160 


27 




2099- WD3 


18300 


900 


8.230 


0.210 


28 




2099- WD4 


16900 


1100 


8.400 


0.260 


29 




2099- WD5 


18300 


1000 


8.330 


0.220 


30 




2099- WD7 


17800 


1400 


8.420 


0.320 


31 
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Table 1 — Continued 



System 



WD name 



Teff [K] a(Teff) \ogig) [cm s"^] a{\og{g)) ID 



2099- WD9 


15300 


400 


8.000 


0.080 


32 


2099- WD 10 


19300 


400 


8.200 


0.070 


33 


2099- WD 11 


23000 


600 


8.540 


0.100 


34 


2099-WD12 


13300 


1000 


7.910 


0.120 


35 


2099- WD 13 


18200 


400 


8.270 


0.080 


36 


2099- WD 14 


11400 


200 


7.730 


0.160 


37 


2099-WD16 


13100 


500 


8.340 


0.100 


38 


NGC2168 LAWDSl 


32400 


512 


8.400 


0.125 


39 


NGC2168 LAWDS 2 


32700 


603 


8.340 


0.080 


40 


NGC2168 LAWDS 5 


52600 


1160 


8.240 


0.095 


41 


NGC2168 LAWDS 6 


55200 


897 


8.280 


0.065 


42 


NGC2168 LAWDS 15 


29900 


318 


8.480 


0.060 


43 


NGC2168 LAWDS 27 


30500 


397 


8.520 


0.061 


44 



NGC2168 (M35) 



Sirius 


Sirius B 


25193 


37 


8.566 


0.010 


45 


NGC7789 


NGC 7789-5 


31213 


238 


7.904 


0.054 


46 




NGC 7789-8 


24319 


447 


8.004 


0.066 


47 




NGC 7789-9 


20939 


727 


7.838 


0.115 


48 


NGC6819 


NGC 6819-6 


21094 


252 


7.832 


0.036 


49 




NGC 6819-7 


15971 


197 


7.908 


0.038 


50 


NGC 1039 


NGC1039 LAWDS 15 


25900 


1100 


8.380 


0.120 


51 




NGC1039 LAWDS 17 


24700 


1100 


8.440 


0.120 


52 




NGC1039 LAWDS S2 


31200 


1100 


8.320 


0.120 


53 
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Table 2. Sources of the data used for determining the cluster ages. 



Name 



[Fe/H] 



E(B-V) 



CMD 



Pleiades 

Hyades 

Praesepe 

NGC2516 

NGC3532 

NGC2099 (M37) 

NGC2168 (M35) 

NGC7789 

NGC6819 

NGC1039 



Gratton (.2000) 
Gratton (2000) 
Gratton (2000 ) 
Gratton f2000) 
Gratton ('2000') 
Chen et al. ('20031 
Gratton ('20001 
Gratton (20001 
Bragadia et al. (20011 
Schuler et al. (20031 



Percival et al. (2003 



Loktin et al. (20011 
Loktin et al. (20011 
Loktin et al. (20011 
Loktin et al. (20011 
Percival fc Salaris (20031 
Bragadia et al. (2001 1 
Loktin et al. (2001) 



Johnson fc Mitchell (19581 
Johnson fc Knuckles (19551 
Johnson (1952) 
Sung et al. (20021 
Fernandez fc Salgado (19801 
Nilakshi fc Sagar (20021 
Sung fc Bessell (19991 
Mocheiska fc Kaluznv (19991 
Kalirai et al. (2001al 
Jones &: Prosser (19961 
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Table 3. Cluster ages determined from the adopted metallicities, reddening values and 
derived distance moduli, using three different sets of isochrones. 



Name 


Fe/H 


E(B- 


-V) 


(m - M)v 


Age[Myr]^ 


Age[Myr]'' 


Age[Myr]'' 




Pleiades 




-0.03 ±0.06 


0.04 




5.74±0.05 


50±10 


85±10 


95±10 


Hyades 




0.13 ±0.06 


0.0 




3.33±0.01 


440±40 


640±40 


630±40 


Praesepe 




0.04 ±0.06 


0.0 




6.24±0.04 


450±40 


650±50 


640±40 


NGC2516 




-0.16 ±0.11 


0.10 




8.27±0.07 


85±45 


130±50 


140±50 


NGC3532 




0.02 ±0.06 


0.04 




8.40±0.25 


300±100 


400±100 


400±100 


NGC2099 


(M37) 'i 


0.09 ±0.15 


0.30 




12.00±0.12 


220±30 


320±30 


320±30 






-0.20 


0.23 




11.40±0.12 


350±40 


550±50 


540±50 


NGC2168 


(M35) 


-0.19 ±0.15 


0.26 




10.50±0.12 


85±25 


120±30 


130±30 


NGC7789 




-0.13 ±0.08 


0.29 




12.12±0.12 


1100±100 


1500±100 


1600±100 


NGC6819 




0.09 ±0.03 


0.14 ±0.04 


12.60±0.20 


1500±200 


2000±200 


2000±200 


NGC1039 




0.07 ±0.04 


0.07 




8.80±0.15 


150±30 


250±25 


250±25 


Sirius A 




0.0 






181±30 


170±50 


170±25 





^Ages from BaSTI isochrones without overshooting 

^Ages from BaSTI isochrones with overshooting 

'^Ages from Padua isochrones with overshooting 

'^For NGC2099 a second alternative [Fe/H] value was tested (see text for details) 
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Table 4. Fractional 1 — a uncertainties in WD mass and cooling age for relevant WD 
input physics, as functions of effective temperature and surface gravity. 





Fract. uncert. for WD mass 


Fract. 
8.6 


uncert 

8.4 


. for WD coolii 

8.2 8.0 


igage 


\og{g) [cm s 2 = 


8.6 


8.4 


8.2 


8.0 


7.8 


7.8 


log(Tcff) [K 
































Coohng 


tracks 










4.60 


0.017 


0.015 


0.016 


0.016 


0.003 


0.826 


0.564 


0.287 


0.335 


0.377 


4.45 


0.016 


0.015 


0.014 


0.010 


0.011 


0.036 


0.262 


0.291 


0.174 


0.006 


4.30 


0.015 


0.014 


0.013 


0.007 


0.017 


0.072 


0.064 


0.012 


0.057 


0.159 


4.15 


0.012 


0.012 


0.011 


0.003 


0.037 


0.035 


0.114 


0.080 


0.029 


0.100 


4.00 


0.009 


0.006 


0.010 


0.000 


0.061 


0.368 


0.020 


0.073 


0.064 


0.122 












Neutrino 


cooling 










4.60 


0.000 


0.001 


0.001 


0.002 


0.003 


0.700 


0.396 


0.320 


0.269 


0.217 


4.45 


0.000 


0.000 


0.000 


0.001 


0.001 


0.147 


0.349 


0.474 


0.374 


0.305 


4.30 


0.000 


0.000 


0.000 


0.000 


0.000 


0.042 


0.075 


0.137 


0.250 


0.326 


4.15 


0.000 


0.000 


0.000 


0.000 


0.000 


0.017 


0.027 


0.044 


0.072 


0.098 


4.00 


0.000 


0.000 


0.000 


0.000 


0.000 


0.007 


0.011 


0.018 


0.027 


0.037 










( 


Ilonductive opacity 








4.60 


0.000 


0.000 


0.000 


0.001 


0.001 


0.391 


0.095 


0.040 


0.022 


0.020 


4.45 


0.000 


0.000 


0.001 


0.001 


0.002 


0.056 


0.197 


0.200 


0.089 


0.046 


4.30 


0.000 


0.000 


0.000 


0.001 


0.002 


0.033 


0.009 


0.028 


0.083 


0.090 


4.15 


0.000 


0.000 


0.000 


0.000 


0.001 


0.056 


0.049 


0.040 


0.024 


0.035 


4.00 


0.000 


0.000 


0.000 


0.000 


0.000 


0.067 


0.058 


0.054 


0.053 


0.083 












Core composition 










4.60 


0.002 


0.002 


0.002 


0.002 


0.002 


0.058 


0.046 


0.043 


0.046 


0.066 


4.45 


0.002 


0.002 


0.002 


0.002 


0.002 


0.013 


0.042 


0.040 


0.043 


0.047 


4.30 


0.002 


0.001 


0.001 


0.002 


0.002 


0.003 


0.028 


0.036 


0.040 


0.043 


4.15 


0.002 


0.001 


0.001 


0.001 


0.002 


0.005 


0.013 


0.033 


0.037 


0.039 


4.00 


0.002 


0.001 


0.001 


0.001 


0.002 


0.001 


0.007 


0.029 


0.037 


0.039 










H-envelope thickness 








4.60 


0.006 


0.008 


0.012 


0.014 


0.020 


0.035 


0.003 


0.011 


0.001 


0.024 


4.45 


0.005 


0.007 


0.010 


0.013 


0.021 


0.029 


0.015 


0.007 


0.004 


0.014 


4.30 


0.005 


0.006 


0.009 


0.012 


0.024 


0.033 


0.011 


0.005 


0.013 


0.006 


4.15 


0.004 


0.006 


0.008 


0.011 


0.027 


0.028 


0.027 


0.007 


0.008 


0.015 
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Table 4 — Continued 



Fract. uncert. for WD mass Fract. uncert. for WD cooling age 



log(^) [cm s-2] = 8.6 8.4 8.2 8.0 7.8 8.6 8.4 8.2 8.0 7.^ 



log(Teff) [K] 



4.00 0.003 0.005 0.007 0.010 0.028 0.070 0.057 0.044 0.025 0.016 

'^Columns and rows are labeled according to log (g) and log (Teg) respectively. Definitions of how 
the uncertainties were computed are given in Eqs. [3H71 Note, however, that for the sake of readability, 
uncertainties given here are linear (not logarithmic) in mass and cooling age and have been assumed 
symmetric. Fractional uncertainties smaller than 0.001 are shown as null values. 
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Table 5. White dwarf cooling ages and associated mass (in solar mass units) determined 
with the SOO cooling tracks. The bolometric luminosity is obtained using the derived 

cluster distances. 



ID 


logtcooi [yr 


o-_(logtcooi) 


o-+(logtcooi) 


Mf[M0] 


a_(Mf) 


^+(Mf) 


log(L/Lo) 




Pleiades 














1 


7.667 


-0.185 


0.162 


1.028 


-0.031 


0.031 


-1.162 


Hyades 














2 


8.262 


-0.063 


0.056 


0.713 


-0.031 


0.031 


-2.032 


3 


8.466 


-0.066 


0.064 


0.798 


-0.032 


0.033 


-2.282 


4 


7.981 


-0.105 


0.086 


0.679 


-0.031 


0.031 


-1.701 


5 


7.585 


-0.185 


0.174 


0.699 


-0.031 


0.031 


-1.324 


6 


7.783 


-0.146 


0.130 


0.652 


-0.032 


0.032 


-1.518 


7 


7.274 


-0.186 


0.184 


0.684 


-0.032 


0.031 


-1.094 


8 


8.425 


-0.063 


0.061 


0.773 


-0.032 


0.033 


-2.237 


Praesepe 














9 


8.113 


-0.072 


0.065 


0.901 


-0.034 


0.034 


-1.739 


10 


8.159 


-0.076 


0.066 


0.622 


-0.028 


0.032 


-1.942 


11 


8.584 


-0.086 


0.085 


0.820 


-0.040 


0.039 


-2.443 


12 


8.338 


-0.064 


0.062 


0.812 


-0.034 


0.034 


-2.087 


13 


8.408 


-0.058 


0.057 


0.857 


-0.023 


0.022 


-2.161 


14 


8.459 


-0.064 


0.063 


0.759 


-0.032 


0.033 


-2.294 


15 


8.498 


-0.062 


0.062 


0.751 


-0.033 


0.033 


-2.357 


16 


8.402 


-0.053 


0.051 


0.722 


-0.026 


0.027 


-2.226 


17 


8.439 


-0.052 


0.051 


0.740 


-0.021 


0.021 


-2.273 


18 


8.476 


-0.061 


0.061 


0.759 


-0.027 


0.027 


-2.320 


19 


8.467 


-0.054 


0.053 


0.746 


-0.027 


0.027 


-2.311 


NGC2516 














20 


7.746 


-0.298 


0.192 


0.926 


-0.106 


0.105 


-1.323 


21 


7.543 


-0.364 


0.297 


1.008 


-0.071 


0.070 


-1.069 


22 


7.888 


-0.116 


0.087 


0.970 


-0.045 


0.045 


-1.455 


23 


7.858 


-0.202 


0.154 


1.075 


-0.078 


0.078 


-1.326 


NGC3532 














24 


7.373 


-0.216 


0.219 


0.482 


-0.085 


0.076 


-1.164 


25 


7.064 


-0.303 


0.161 


0.563 


-0.169 


0.145 


-0.789 


26 


8.047 


-0.284 


0.213 


0.708 


-0.150 


0.164 


-1.762 
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Table 5 — Continued 



ID logtcooi [yr] a_(logteooi) a+(logteooi) MffMo] a_(Mf) a+(Mf) \og{L/Lrs 



NGC2099 (M37) 














27 


7.971 


-0.215 


0.150 


0.690 


-0.093 


0.096 


-1.685 


28 


8.185 


-0.193 


0.167 


0.759 


-0.128 


0.133 


-1.909 


29 


8.411 


-0.214 


0.245 


0.864 


-0.163 


0.164 


-2.161 


30 


8.260 


-0.190 


0.194 


0.821 


-0.135 


0.139 


-1.975 


31 


8.366 


-0.268 


0.305 


0.877 


-0.198 


0.203 


-2.085 


32 


8.264 


-0.083 


0.072 


0.614 


-0.044 


0.051 


-2.082 


33 


8.090 


-0.094 


0.074 


0.742 


-0.044 


0.044 


-1.796 


34 


8.111 


-0.108 


0.097 


0.959 


-0.065 


0.064 


-1.720 


35 


8.380 


-0.133 


0.120 


0.561 


-0.080 


0.067 


-2.275 


36 


8.222 


-0.082 


0.075 


0.783 


-0.050 


0.052 


-1.945 


37 


8.455 


-0.125 


0.121 


0.448 


-0.129 


0.123 


-2.461 


38 


8.667 


-0.097 


0.103 


0.818 


-0.064 


0.066 


-2.567 


NGC2168 (M35) 














39 


7.343 


-0.339 


0.355 


0.884 


-0.078 


0.079 


-1.021 


40 


7.223 


-0.310 


0.327 


0.849 


-0.051 


0.050 


-0.962 


41 


6.351 


-0.213 


0.208 


0.820 


-0.055 


0.056 


-0.052 


42 


6.292 


-0.218 


0.217 


0.845 


-0.041 


0.040 


+0.006 


43 


7.634 


-0.219 


0.200 


0.928 


-0.040 


0.040 


-1.219 


44 


7.644 


-0.215 


0.186 


0.955 


-0.041 


0.041 


-1.212 


Sirius 


B 














45 


7.996 


-0.063 


0.062 


0.978 


-0.016 


0.016 


-1.580 


NGC7789 














46 


6.998 


-0.137 


0.135 


0.601 


-0.028 


0.029 


-0.757 


47 


7.480 


-0.181 


0.179 


0.637 


-0.037 


0.039 


-1.265 


48 


7.627 


-0.205 


0.192 


0.543 


-0.059 


0.058 


-1.428 



NGC6819 














49 7.609 


-0.153 


0.150 


0.540 


-0.021 


0.019 


-1.412 


50 8.141 


-0.071 


0.063 


0.566 


-0.023 


0.021 


-1.951 



NGC1039 
51 7.801 



-0.239 



0.161 



0.863 



-0.075 



0.075 



-1.443 
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Table 5 — Continued 



ID 


logtcool M 0-_(logtcool) 


0-+(l0gtcool) 


Mi[Mq 


a_(Mf) 


^+(Mf) 


log(L/Lo) 


52 
53 


7.934 
7.308 


-0.179 
-0.312 


0.126 
0.339 


0.899 
0.834 


-0.076 
-0.072 


0.076 
0.076 


-1.533 

-1.144 
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Table 6. Initial-final-mass- relation obtained by combining BaSTI progenitor models with 
and without overshooting (indicated by OV resp. non-OV) with SOO WD tracks ^ . 





Mf[MQ 


a(Mf) 




OV models 




non 


-OV models 


ID 


Mi [Mo] 


a- (Mi) 


a+(Mi) 


Mi [Mo] 


^-(Mi) 


o-+(Mi) 




Pleiades 
















1 


1.028 


0.031 


7.847 


-1.257 


3.069 


15.52 


-3.986 


99.90 


Hyades 
















2 


0.713 


0.031 


3.067 


-0.100 


0.107 


3.690 


-0.192 


0.224 


3 


0.798 


0.032 


3.383 


-0.174 


0.208 


4.438 


-0.454 


0.813 


4 


0.679 


0.031 


2.887 


-0.076 


0.079 


3.385 


-0.128 


0.142 


5 


0.699 


0.031 


2.771 


-0.067 


0.066 


3.184 


-0.108 


0.111 


6 


0.652 


0.032 


2.815 


-0.071 


0.071 


3.261 


-0.116 


0.122 


7 


0.684 


0.031 


2.746 


-0.063 


0.064 


3.144 


-0.097 


0.105 


8 


0.773 


0.032 


3.284 


-0.151 


0.172 


4.240 


-0.349 


0.549 


Praesepe 
















9 


0.901 


0.034 


2.799 


-0.094 


0.106 


3.292 


-0.138 


0.161 


10 


0.622 


0.030 


2.827 


-0.099 


0.112 


3.341 


-0.149 


0.179 


11 


0.820 


0.039 


3.543 


-0.337 


0.548 


5.953 


-1.223 


8.634 


12 


0.812 


0.034 


2.972 


-0.130 


0.148 


3.655 


-0.230 


0.287 


13 


0.857 


0.023 


3.058 


-0.153 


0.174 


3.926 


-0.290 


0.421 


14 


0.759 


0.032 


3.179 


-0.183 


0.221 


4.151 


-0.400 


0.696 


15 


0.751 


0.034 


3.282 


-0.206 


0.260 


4.457 


-0.524 


1.048 


16 


0.722 


0.027 


3.045 


-0.146 


0.167 


3.901 


-0.272 


0.386 


17 


0.740 


0.021 


3.130 


-0.160 


0.188 


4.060 


-0.328 


0.502 


18 


0.759 


0.027 


3.223 


-0.189 


0.232 


4.233 


-0.447 


0.801 


19 


0.746 


0.027 


3.200 


-0.175 


0.210 


4.190 


-0.397 


0.654 


NGC2516 
















20 


0.926 


0.106 


5.856 


-1.133 


3.679 


8.672 


-2.774 


99.90 


21 


1.008 


0.071 


5.269 


-0.837 


1.899 


6.732 


-1.511 


5.739 


22 


0.970 


0.045 


6.766 


-1.541 


7.174 


13.99 


-5.870 


99.90 


23 


1.075 


0.078 


6.492 


-1.415 


5.640 


12.35 


-4.929 


99.90 


NGC3532 
















24 


0.482 


0.081 


3.132 


-0.252 


0.369 


3.450 


-0.340 


0.563 


25 


0.563 


0.157 


3.085 


-0.243 


0.348 


3.409 


-0.321 


0.523 
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Table 6 — Continued 





MffM©; 


a(Mf) 




OV models 




non 


-OV models 


ID 


M,[Mq] 


a_(MO 


<J+{M,) 


M;[Mq 


a4Mi) 


cr+{M-^ 


26 


0.708 


0.157 


3.466 


-0.399 


0.620 


3.975 


-0.560 


1.316 


NGC2099 (M37) 


i; [Fe/H 


= 0.09 












27 


0.690 


0.095 


3.934 


-0.248 


0.326 


4.746 


-0.491 


0.868 


28 


0.759 


0.131 


4.372 


-0.436 


0.940 


6.009 


-1.110 


6.012 


29 


0.864 


0.164 


6.307 


-1.557 


99.90 


99.90 


99.90 


99.90 


30 


0.821 


0.137 


4.666 


-0.567 


1.891 


7.628 


-2.020 


99.90 


31 


0.877 


0.200 


5.497 


-1.099 


99.90 


99.90 


99.90 


99.90 


32 


0.614 


0.048 


4.689 


-0.437 


0.693 


7.806 


-1.834 


8.624 


33 


0.742 


0.044 


4.153 


-0.238 


0.305 


5.187 


-0.593 


1.015 


34 


0.959 


0.065 


4.198 


-0.274 


0.382 


5.323 


-0.691 


1.433 


35 


0.561 


0.074 


5.680 


-1.022 


4.981 


99.90 


99.90 


99.90 


36 


0.783 


0.051 


4.471 


-0.364 


0.519 


6.574 


-1.206 


3.976 


37 


0.448 


0.126 


8.210 


-2.545 


99.90 


99.90 


99.90 


99.90 


38 


0.818 


0.065 


99.90 


99.90 


99.90 


99.90 


99.90 


99.90 


NGC2168 (M35) 


1 














39 


0.884 


0.079 


5.209 


-0.595 


0.928 


6.119 


-0.915 


1.695 


40 


0.849 


0.050 


5.111 


-0.522 


0.796 


5.936 


-0.784 


1.330 


41 


0.820 


0.055 


4.903 


-0.387 


0.605 


5.468 


-0.550 


0.890 


42 


0.845 


0.040 


4.899 


-0.385 


0.602 


5.461 


-0.549 


0.885 


43 


0.928 


0.040 


5.767 


-0.862 


1.670 


7.312 


-1.527 


4.464 


44 


0.955 


0.041 


5.804 


-0.869 


1.669 


7.383 


-1.558 


4.441 


Sirius B 
















45 


0.978 


0.016 


5.934 


-1.139 


3.904 


5.426 


-0.777 


1.198 


NGC7789 
















46 


0.601 


0.029 


1.853 


-0.096 


0.128 


1.954 


-0.072 


0.105 


47 


0.637 


0.038 


1.864 


-0.100 


0.132 


1.971 


-0.078 


0.116 


48 


0.543 


0.059 


1.871 


-0.101 


0.134 


1.983 


-0.083 


0.125 


NGC6819 
















49 


0.540 


0.020 


1.736 


-0.090 


0.117 


1.852 


-0.073 


0.077 


50 


0.566 


0.022 


1.770 


-0.097 


0.130 


1.887 


-0.072 


0.096 
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Table 6 — Continued 





Mf[MQ 


a(Mf) 




OV models 




non 


-OV models 


ID 


m^iMq] 


^-(Mi) 


cr^M,) 


M;[Mq] 


a_{MO 


^+iM,) 




NGC1039 
















51 


0.863 


0.075 


4.124 


-0.364 


0.514 


5.382 


-0.717 


1.369 


52 


0.899 


0.076 


4.319 


-0.435 


0.659 


6.069 


-1.006 


2.652 


53 


0.834 


0.074 


3.835 


-0.279 


0.348 


4.607 


-0.386 


0.511 


NGC2099 (M37; 


); [Fe/H] = 


-0.20 












27 


0.690 


0.095 


2.994 


-0.133 


0.145 


3.643 


-0.269 


0.349 


28 


0.759 


0.131 


3.123 


-0.196 


0.259 


4.029 


-0.443 


0.947 


29 


0.864 


0.163 


3.522 


-0.355 


1.056 


5.248 


-1.460 


99.90 


30 


0.821 


0.137 


3.232 


-0.233 


0.374 


4.226 


-0.570 


1.727 


31 


0.877 


0.201 


3.424 


-0.325 


1.085 


4.814 


-1.064 


99.90 


32 


0.614 


0.047 


3.239 


-0.173 


0.200 


4.238 


-0.486 


0.844 


33 


0.742 


0.043 


3.055 


-0.128 


0.144 


3.823 


-0.273 


0.382 


34 


0.959 


0.065 


3.068 


-0.139 


0.156 


3.865 


-0.306 


0.443 


35 


0.561 


0.074 


3.454 


-0.282 


0.412 


4.923 


-1.010 


4.628 


36 


0.783 


0.051 


3.174 


-0.162 


0.181 


4.121 


-0.408 


0.645 


37 


0.448 


0.125 


3.628 


-0.361 


0.628 


6.066 


-2.242 


99.90 


38 


0.818 


0.066 


5.529 


-1.125 


9.156 


99.90 


99.90 


99.90 



^Values for initial masses equal to 99.90 M0 flag those cases where the central value of the 
progenitor age tprog, is negative. Similarly, error values of 99.90 indicate that progenitor ages 
deviating l-a from the central value tprog are negative and as a consequence l-a values for Mi 
cannot be obtained. 
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Fig. 1. — Comparison of the cluster ages determin ed with the Ba S TI ov ershooting isochrones, 



(120051), for the clusters in 



with the ages adopted from the literature by iFerrario et al. 

common. The dash ed line displays the 1:1 relation, and the errors are those given by 



Ferrario et al.l (120051 ) or as listed from Table [31 The open circle shows the results for Sirius A. 
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Fig. 2. — Reference IFMR obtained employing BaSTI isochrones and progenitor models with 
core overshooting together with SCO WD tracks. Lines show fits to the data: the solid line is 
a linear fit, the dashed line a piecewise linear fit. The pivot point at 4 Mq has been chosen 
following theoretical predictions from BaSTI models, that show changes in the slope of the 
IFMR around those mass values. 
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Fig. 3. — Inferred and theoretical IFMRs. Top panel: results for models and isochrones 
with core overshooting. Solid lines correspond to theoretical IFMRs derived from BaSTI 
models (thick lines - for [Fe/H]=0.06) and LPCODE models (thin line - for [Fe/H]= 0.06). 
The dotted thick line shows the M; — MdTp relation for BaSTI models and the dashed thin 
line shows the counterpart from LPCODE models. Bottom panel: results corresponding to 
models and isochrones without core overshooting. In this case only theoretical results for 
BaSTI models are shown. 
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Fig. 4. — Difference between initial masses obtained employing the Padua and BaSTI 
isochrones and models versus initial masses obtained with BaSTI isochrones and models. 
Both sets include overshooting treated with different formulations, but giving similar mass 
extensions of the core overshooting regions at fixed total stellar mass. SOO WD tracks have 
been used in both cases. 
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Fig. 5. — Comparison between IFMRs obtained using the SOO (dots) and the LPCODE 
(triangles) WD tracks. BaSTI isochrones and progenitor models with core overshooting 
have been used in both cases. Each pair connected by a line represents the same star. 
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Fig. 6. — Contribution of uncertainty sources to the total error budget for WD mass (upper 
left), WD age (upper right), progenitor mass (lower left), and progenitor age (lower right 
panel). The fractional contribution to total sigma a due to various error sources is given 
for each WD. Vertical lines separate clusters. For the WDs we show the importance of 
observational vs. input physics and systematics (different codes, variation of physics inputs, 
chemical stratification, envelope thickness) for the progenitors that of cluster and WD age. 
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Fig. 7.— IFMR for NGC 2099, for the two [Fe/H] and E(B-V) pairs given in Table El Filled 
circles refer to [Fe/H]= 0.09 and empty circles to [Fe/H]= —0.20. The solid line shows the 
theoretical IFMR from BaSTI models with core overshooting, the dotted line displays the 
Mi — MciTp relation for the same models. Data points lie systematically below the theoretical 
Mi — MciTp values when [Fe/H]= 0.09 is adopted. 



